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Abstract 

As  an  emerging  effective  approach  to  nonlinear  robust  control,  simplex  sliding  mode  control  demonstrates  some 
attractive  features  not  possessed  by  the  conventional  sliding  mode  control  method,  from  both  theoretical  and  practical 
points  of  view.  However,  no  systematic  approach  is  currently  available  for  computing  the  simplex  control  vectors  in 
nonlinear  sliding  mode  control.  In  this  paper,  chaos-based  optimization  is  exploited  so  as  to  develop  a  systematic 
approach  to  seeking  the  simplex  control  vectors;  particularly,  the  flexibility  of  simplex  control  is  enhanced  by  making 
the  simplex  control  vectors  dependent  on  the  Euclidean  norm  of  the  sliding  vector  rather  than  being  constant,  which 
result  in  both  reduction  of  the  chattering  and  speedup  of  the  convergence.  Computer  simulation  on  a  nonlinear 
uncertain  system  is  given  to  illustrate  the  effectiveness  of  the  proposed  control  method. 

©  2004  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

As  a  powerful  tool  for  solving  complex  control  problems  under  significant  uncertainties,  sliding  mode  control 
(SMC)  has  been  widely  studied  and  applied  in  different  fields  of  engineering  applications  over  the  last  30  years.  The 
primary  feature  of  sliding-mode  control  systems  is  that  the  feedback  signal  is  discontinuous,  switching  on  one  or 
more  manifolds  in  the  state  space.  When  the  system  state  crosses  each  discontinuity  surface,  the  structure  of  the 
feedback  system  is  altered.  Under  certain  circumstances,  all  motions  in  the  neighborhood  of  the  desired  manifold  are 
directed  towards  the  manifold  and  thus  a  sliding  motion  on  a  predefined  subspace  within  the  state  space  is  estab¬ 
lished,  in  which  the  system  state  repeatedly  crosses  the  switching  surface.  During  the  sliding  mode,  the  system 
possesses  high  robustness  against  uncertainties  of  various  kinds,  such  as  parameter  variations  and  external  distur¬ 
bances. 

However,  due  to  the  use  of  discontinuous  control  signals,  the  sliding  mode  control  inherently  suffers  from  the 
chattering  problem.  For  an  m-input  system,  the  conventional  sliding-mode  control  approach  usually  partitions  the  state 
space  into  2m  regions  and  the  control  law  will  switch  whenever  the  system  state  moves  from  one  region  to  another.  As  a 
consequence,  unwanted  high-frequency  chattering  motions  are  generated,  which  degrades  the  system  performance  in 
general.  To  reduce  the  switching  behaviour  and  circumvent  the  redundancy  of  control  actions  that  occurs  in  the  case  of 
the  conventional  component- wise  sliding-mode  control  design  procedure,  Baida  and  Izosimov  [1]  proposed  a  new 
algorithm,  called  the  simplex  sliding  mode  control  (SSMC)  scheme,  which  only  partitions  the  state  space  into  m  +  1 
regions.  Clearly,  the  number  of  regions  is  decreased  and  as  a  result  the  chattering  problem  can  be  improved.  Besides, 
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SSMC  is  also  an  effective  way  to  extend  the  sliding-mode  control  methodology  to  the  multi-input  case.  The  design  of  the 
so-called  simplex  control  for  multivariable  systems  [2,3]  is  relatively  straightforward  in  comparison  with  the  conven¬ 
tional  SMC,  and  some  successful  applications  of  SSMC  in  practice  have  been  reported  [4]. 

SSMC  differs  in  certain  fundamental  aspects  from  the  conventional  SMC.  For  instance,  the  state  space  is  partitioned 
into  only  m  +  1  regions,  instead  of  2m  regions  as  is  typical  for  conventional  SMC,  where  m  is  the  number  of  control 
inputs.  Furthermore,  a  distinct  control  vector  is  used  in  each  region  of  the  partition  for  the  simplex  control  scheme,  and 
these  m  +  1  control  vectors  satisfy  the  simplex  property  (to  be  defined  later)  in  a  reduced-order  subspace  whose  origin 
corresponds  to  the  intersection  of  the  surfaces  Spc  =  0,  i  =  1 , . . . ,  m,  which  is  the  desired  sliding  mode  hypersurface  in 
the  state  space.  In  addition,  for  this  scheme  (except  for  m  =  1)  the  surfaces  of  control  discontinuity  do  not  coincide  with 
the  linear  sliding  hypersurfaces  Stx  =  0,  i  =  1 , . . . ,  m,  as  is  the  case  for  conventional  SMC;  instead,  there  exist 
m(m  +  l)/2  switching  hypersurfaces  for  the  control. 

In  other  respects,  the  SSMC  design  methodology  is  similar,  in  principle,  to  the  conventional  SMC  design.  Having 
first  selected  S  =  [Sff  Sj  •  •  •  S^]T  to  obtain  some  desired  sliding-mode  dynamics  on  the  hypersurface  Sx  =  0 
(sliding  phase  design),  choose  next  the  control  vectors  so  that  the  closed-loop  system  trajectory  will  be  forced  into  a 
sliding  mode  on  that  hypersurface  (reaching  phase  design). 

The  simplicity  and  applicability  of  SSMC  are  so  appealing  that  there  have  been  some  publications  about  how  to 
synthesize  the  SSMC  law  for  linear  time-invariant  systems  [5,6];  however,  no  systematic  approach  is  currently  available 
for  computing  the  simplex  control  vectors  in  nonlinear  sliding  mode  control.  In  this  paper,  an  innovative  global 
optimization  method,  chaos  optimization,  is  exploited  to  develop  a  systematic  approach  to  seeking  the  simplex  control 
vectors  for  nonlinear  sliding  mode  control  systems,  which  significantly  extends  the  state-of-the-art  advances  of  the 
methodology. 

The  rest  of  the  paper  is  organized  as  follows.  In  the  next  section,  the  simplex  sliding  mode  control  and  the  chaos 
optimization  are  introduced.  Section  3  presents  the  procedure  for  synthesizing  the  simplex  sliding  mode  control  law  via 
the  chaos  optimization.  In  Section  4,  the  proposed  synthesis  approach  is  simulated  on  an  uncertain  multi-input  non¬ 
linear  system  to  confirm  its  validity.  Finally,  some  concluding  remarks  are  given  in  Section  5. 


2.  Preliminaries 

2.1.  Simplex-type  nonlinear  sliding-mode  control 
Considering  a  class  of  affine  nonlinear  systems: 

X  =  f(x)  +  g(x)u,  (1) 

where  x  e  Rn,u  e  Rm.  Choose  a  manifold  described  by  s(v)  =  Hx  =  0,  where  s(x)  e  Rm,  so  that  if  the  system  trajectories 
lie  on  the  manifold,  the  behaviour  of  the  system  satisfies  a  prespecified  control  objective.  Assume  that  a  control  signal  is 
designed  so  that  it  is  capable  of  forcing  the  system  trajectories  from  any  initial  state  to  reach  the  manifold  in  finite  time 
and  then  to  slide  on  it. 

Sliding-mode  control  (SMC)  in  the  single-input  case  is  quite  straightforward  and  has  been  well  developed.  In 
contrast,  most  methods  of  designing  multi-input  SMC  for  nonlinear  uncertain  systems  are  far  from  a  simple  gener¬ 
alization  of  the  standard  single-input  variable  structure  approach. 

In  order  to  provide  a  control  strategy  accomplishing  the  objective  of  steering  s(x)  to  zero  in  the  multi-input  case, 
consider  the  dynamic  behaviour  of  the  m-vector  s(x)  associated  with  system  (1),  namely, 

s(x)  =  Hf(x)  +  Hg(x)u ,  (2) 

where  the  matrix  [//g(x)]  is  assumed  to  be  always  invertible. 

A  traditional  way  of  extending  the  basic  SMC  strategy  to  the  multi-input  case  consists  of  selecting  the  control  signals 
according  to  a  component- wise  procedure.  This  means  that  the  control  is  designed  so  that  each  component  of  the  vector 
s(x)  is  related  to  the  corresponding  time  derivative  s(x)  by  the  so-called  reaching  condition.  Consider  for  instance  the 
constant  plus  proportional  rate  reaching  law  [7]: 

s  =  —  2  sign (5)  -Ks,  (3) 

where  Q  =  diagjgf },  K  =  diag{^},  i  =  1, . . . ,  m  function  as  the  weighting  matrices.  The  control  law  derived  from  (2) 
and  (3)  becomes 
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u  =  -  [Hg{x)  r1  [Hf(x)  +  gsignCs)  +  Ks] .  (4) 

It  is  noteworthy  that  each  component  of  the  vector  u  is  a  linear  combination  of  some  discontinuous  signals. 

An  alternative  approach  to  multi-input  nonlinear  sliding-mode  control  is  the  simplex  control  [2,3].  The  advantage  of 
the  simplex  control  approach  lies  in  the  fact  that  the  number  of  structures  among  which  the  controlled  system  switches 
is  reduced  with  respect  to  the  component-wise  case  (m  +  1  instead  of  2m).  To  introduce  this  method,  the  following 
definitions  are  needed. 


Definition  1.  The  set  U  =  {V,  w2, . . . ,  um+l},  where  ul  e  Rm,  i  =  1, 2, . . . ,  m  +  1,  are  distinct  and  nonzero  vectors,  is  said 
to  form  a  simplex  in  Rm  if  every  subset  of  m  vectors  in  U  are  linearly  independent  and  there  exist  m  +  1  real  and  positive 
constants  {ai,  a2, . . . ,  am+i}  such  that 

m+ 1  m+1 

a.  =  i  and  a  tul  =  0.  (5) 

*=i  »=i 


This  definition  means  that  a  simplex  is  a  set  of  m  +  1  affinely  independent  vectors  in  Rm  such  that  0m  is  in  the  interior  of 
the  convex  hull  of  those  vectors.  For  simplex-type  nonlinear  sliding-mode  control,  the  following  transformation  is 
performed: 

<7=  [/fg(x)]_1^(x).  (6) 

Next,  define  the  m+1  open  regions  in  the  m- dimensional  cr-space  associated  with  such  a  simplex  as 

=  |  a  :  <7  =s  2jUJ\  Xj  >  0 1,  i  =  1 , 2, . . . ,  m  +  1 .  (7) 

l  j—  i  j+  J 

Thus,  each  is  an  infinite  cone,  with  vertex  at  0m,  situated  on  the  side  of  Rm  opposite  to  u\ 

Definition  2.  Suppose  that  the  set  U  =  {ul,  u2, . . .  ,wm+1}  forms  a  simplex  in  the  m-dimensional  cr-space.  The  simplex 
control  law  is  defined  to  be 


u  =  u(a)  =  u\  for  o  G  (8) 

Note  that  given  a  simplex  U,  the  control  switching  surfaces  for  the  simplex  control  law  are  defined  by  the  boundaries  of 
Qh  which  is  apparently  not  accordant  with  the  sliding  manifold.  The  following  theorem  [2,3]  investigates  the  stability 
and  convergence  of  the  simplex  nonlinear  sliding  mode  control  approach. 

Theorem  1.  Partition  the  o-state  space  into  m  +  1  non-overlapping  regions  as  (7).  If  the  control  is  defined  by  (8),  then  the 
origin  o  =  0  is  asymptotically  stable  and  the  convergence  takes  place  in  finite  time. 

From  (6),  it  is  clear  that  the  above  theorem  guarantees  the  closed-loop  system  defined  by  (1)  and  (8)  will  achieve 
the  sliding  mode  in  finite  time,  i.e.,  s(x)  =  Hx  =  0.  Note  that  the  control  is  unique  in  each  region,  and  each  region  is 
uniquely  determined  by  the  choice  of  the  simplex  of  control  vectors. 

2.2.  Chaos  optimization 

Over  the  last  decade,  one  has  been  experienced  a  resurgence  of  interest  in  the  control  systems  society  for  new 
methods  of  global  optimization  as  well  as  the  application  of  available  global  optimization  algorithm  to  control  systems 
analysis  and  synthesis.  This  recent  surge  of  interest  is  attributed  to  the  fact  that  a  number  of  problems  in  the  field  of 
robust  control  may  be  directly  restated  as  an  optimization  problem. 

Genetic  algorithm  (GA)  and  simulated  annealing  (SA)  are  among  the  most  popular  tools  for  global  optimization 
[8,9].  By  mimicking  the  metaphor  of  natural  biological  evolution,  GA  operates  on  a  population  of  potential  solutions 
applying  the  principle  of  survival  of  the  fittest  to  produce  successively  better  approximations  to  a  solution.  Inspired  by 
the  annealing  process  for  metals  during  cooling,  SA  describes  a  family  of  iterative  methods  where  every  iteration 
consists  of  taking  a  step  in  the  parameter  space  with  a  probability  that  decays  exponentially  with  an  “energy”  function 
associated  with  the  new  parameter  value. 
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As  derivative-free  stochastic  search  methods  for  determining  the  solutions  of  an  optimization  problem,  simulated 
annealing  and  evolutionary  computation  both  employ  a  stochastic  mechanism  to  avoid  being  trapped  in  a  local 
optimum.  On  the  contrary,  the  deterministic  dynamics  plays  a  central  role  in  the  method  of  chaos  optimization  [10]. 

Chaos,  an  apparently  disordered  behavior  that  is  nonetheless  deterministic,  is  a  universal  phenomenon  that  occurs  in 
many  nonlinear  systems.  It  is  featured  by  highly  unstable  motion  of  deterministic  systems  in  a  bounded  region  of  the 
phase  space.  High  instability  means  that  the  distance  of  two  nearby  orbits  increases  exponentially  with  time  [11],  which 
is  a  result  of  the  extreme  sensitivity  of  chaotic  systems  to  the  initial  conditions.  The  Lyapunov  exponents  quantify  this 
property.  The  magnitude  of  a  Lyapunov  exponent  represents  the  principal  rate  of  orbit  divergence  in  the  phase  space. 
For  a  one-dimensional  dynamical  system  xi+\  =  /(v*),  the  Lyapunov  exponent  X  is  defined  as  the  long-time  average  of 
the  exponent  with  respect  to  an  orbit: 


A  =  iim  n  23  los  i /'(*<•) 

jV— X50  TV  z ' 


(9) 


Chaos  is  then  characterized  by  the  boundedness  of  the  system  trajectories  with  a  positive  Lyapunov  exponent,  which 
implies  that  the  average  gradient  of  the  map  is  greater  than  unity,  and  accordingly  two  nearby  orbits  in  phase  space 
diverge  at  an  exponential  rate. 

It  was  emphasized  that  the  sensitivity  to  the  initial  value  suggests  the  irregularity  of  the  series  {v,}  generated  by  chaos 
[11].  Consider  the  ith  number  xt  of  the  series.  It  may  be  possible  that  Xj  with  j  >  i  is  quite  close  to  xf.  Unless  xt  =  Xj 
exactly,  however,  the  part  xuxi+ \,xi+2, ...  is  very  different  from  the  part  Xj,xj+\,xj+2, . . .  due  to  the  sensitivity  to  the 
initial  difference. 

Although  the  long-term  behaviour  of  a  chaotic  system  shows  typical  stochastic  properties,  chaos  is  not  equivalent  to 
a  random  process.  A  chaotic  motion  can  traverse  every  state  in  a  certain  region  (called  the  chaos  space)  by  its  own 
regularity,  and  every  state  is  visited  only  once  therefore  no  precise  periodicity.  The  unique  ergodicity  and  the  irregu¬ 
larity  of  the  series  generated  by  chaos  make  chaotic  dynamics  a  potential  candidate  in  the  field  of  global  optimization 
[10].  In  fact,  it  has  been  successfully  applied  in  solving  nonlinear  programming  problems  [12]  and  improving  the 
performance  of  the  GA  [13]. 

In  our  study,  the  logistic  mapping  is  used  in  chaos  optimization  to  generate  the  chaotic  time  series.  Consider  the 
equation  of  logistic  mapping 

*/+i  =  4xf(l  -xi),  (10) 

where  xt  is  the  chaotic  variable.  Its  Lyapunov  exponent  is 

X  =  log  2  >  0.  (11) 

It  is  well  known  that  chaotic  evolutions  could  be  generated  by  Eq.  (10),  and  the  ergodic  area  (i.e.,  chaos  space)  is  the 
interval  (0, 1).  A  general  procedure  of  chaos  optimization  can  be  found  in  [10]. 


3.  Synthesis  of  the  simplex  control  law  via  chaos  optimization 

The  simplex  sliding  mode  control  is  based  on  a  set  of  m  +  1  control  vectors  forming  a  simplex  in  Rm,  and  on  the 
corresponding  switching  of  the  controlled  system  from  one  to  another  among  m  +  1  different  structures.  In  other  words, 
this  approach  consists  of  the  selection  of  a  set  of  control  vectors  forming  a  simplex,  along  with  a  suitable  switching 
logic.  Hence,  the  crux  of  simplex  control  is  the  choice  of  a  proper  set  of  simplex  control  vectors.  Once  they  are 
appropriately  decided,  the  system  state  space  will  be  partitioned  into  m  +  1  non-overlapping  regions,  and  with  each 
region  a  particular  control  vector,  among  those  of  the  simplex,  is  associated  in  such  a  way  that  the  system  trajectory  is 
forced  to  slide  on  a  prespecified  manifold. 

In  our  study,  the  simplex  control  vectors  were  sought  via  chaos  optimization.  Contrary  to  the  use  of  constant 
simplex  in  the  previous  works,  the  simplex  control  vectors  employed  here  is  dependent  on  the  Euclidean  norm  of  the 
sliding  vector  so  that  the  control  effort  is  proportional  to  the  energy  of  the  sliding  vector,  which  results  in  reduction  of 
the  chattering  and  speedup  of  the  convergence.  Clearly,  the  simplex  control  vectors  include  a  total  of  m  x  (m  +  1) 
entries,  where  m  is  the  number  of  system  inputs. 

It  turns  out  to  be  much  easier  to  describe  the  proposed  method  by  working  out  a  concrete  example  in  detail,  rather 
than  using  some  general  statements  with  complicated  formulas.  For  this  purpose,  consider  the  following  uncertain 
nonlinear  system: 
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x\  =  a{x2  —  Jti), 

x2  =  (c  —  a)x  i  —  x\x2  +  cx2 , 

X3  =  X\X2  —  bx  3, 
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where  a  =  45,  b  E  [3.5, 11],  c  E  [25,  31].  This  nonlinear  model  is  based  on  the  chaotic  Chen’s  system  [14],  and  it  follows 
from  its  bifurcation  analysis  [15]  that  this  interval  system  remains  to  be  chaotic  when  c  is  fixed  to  be  28.  It  has  been 
widely  experienced  that  this  chaotic  system  is  relatively  difficult  to  control  as  compared  to  the  Lorenz  system  and 
Chua’s  system,  which  are  all  topologically  non-equivalent,  due  to  the  prominent  three-dimensional  and  complex 
topological  features  of  its  attractor,  especially  its  rapid  change  in  velocity  in  the  x3 -direction. 

The  interval  system  (12)  can  be  represented  by  a  simple  state-space  equation  of  the  form 

x(t)  =fo(x)  +  A f(x),  x(t)  G  R\  (13) 

where  i;  =  f0(x)  the  nominal  system,  which  can  be  described  by 

{*i  a0(x2  -xi), 

*2  =  (c0  -  a0)xt  -  xxx3  +  c0x2,  (14) 

=  X1X2  -  Z?0X3, 


where  a0  =  45,  Z?0  =  7.25,  c0  =  28  and  A f(x)  denotes  the  system  perturbations.  The  controlled  uncertain  nonlinear 
system  is 

x(t)  =  f0(x)  +  A f{x)  +Bu(t),  (15) 


where  u(t)  e  R2  and 


Given  b  =  bo  +  Ab  and  c  =  c0  +  Ac,  the  uncertainties  can  be  decomposed  into 

a/w.  = f i »)  ■ 


(16) 


Obviously,  the  uncertainties  A f(x)  in  system  (15)  satisfy  the  matching  conditions ,  namely,  A f(x)  is  in  the  subspace 
spanned  by  the  column  vectors  of  B.  So,  the  system  in  sliding  mode  is  robust  to  the  system  perturbations  A f(x). 

In  order  to  find  a  sliding  mode  control  law,  u(t)  E  R2,  which  can  guide  the  system  state  x(t)  to  track  the  prespecified 
reference  signal  xr(t),  a  sliding  surface  needs  to  be  defined  for  the  nominal  system  beforehand.  The  desired  sliding 
manifold  can  be  obtained  by  using  the  method  proposed  in  [16]: 

S(e)  =  {e\s(ejt)  =  He(t)  =  0},  s(e)  E  i?2,  (17) 

where  e  =  x  —  xr  is  the  tracking  error  vector  for  the  desired  trajectory  xr(t). 

Then,  one  can  start  the  procedure  for  synthesizing  the  simplex  control  vectors.  Firstly,  the  performance  index  was 
defined  as 


nu)  =  yyiMi+jU||/M},  (is) 

i=i 

where  N  is  the  duration  of  the  simulation  for  evaluating  the  design,  i  the  time  index  in  simulation,  sf  the  sliding  vector  at 
simulation  step  i,  A st  the  change  of  the  sliding  vector  and  a  bias  weighting  factor  between  st  and  A sh  The  A st  term  can 
be  distinctively  weighted  to  further  suppress  oscillations.  The  simplex  control  vectors  will  be  decided  to  minimize  the 
performance  index  (18),  where  ||  •  ||  denotes  the  Euclidean  norm  of  a  vector. 

Since  six  components  are  involved  in  the  simplex  control  vectors  U  =  { ul,u2,u3},  six  initial  chaotic  variables, 
7i  o,  y2;o,  •  •  • ,  76,o?  °<77-o<  1,7  =  1, 2, . . . ,  6,  are  selected  randomly,  and  the  fixed  points  0.25,  0.5,  0.75  of  the  logistic  map 
cannot  be  used  as  initial  variables.  The  lower  bounds  and  upper  bounds  of  the  searched  variables  are  denoted  as  lbd7 
and  ubd7,  j  =  1,2, . . . ,  6,  and  U*  and  T*  are  assumed  to  be  the  optima.  The  search  procedure  by  chaos  optimization 
starts  from  n  =  0: 
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Step  1:  Chaotify  the  variables: 

Substitute  yX  n,  y2,w,  •  •  • ,  y6,«  the  equation  of  (10)  to  generate  six  chaotic  variables  yln+1,  72,«+o  •  •  •  >  76,«+i  yia  the 
logistic  map. 

Ste/?  2:  Perform  the  transformation  from  the  chaotic  space  to  the  solution  space  by  using  the  following  formula: 

Pj,n+ 1  =  Ibd,  +  ( Libel,  -  lbd,)  X  W,  j  =  1, 2, ...  ,6.  (19) 


Step  3:  Determine  if  K  =  {/c1,  k2,  k3},  where  k1  =  [ftX2-i,«+n  Ax2,«+i]T>  from  a  convex  hull  and  the  origin  of  a-space  is 
enclosed  by  the  convex  hull,  i.e.,  O2  E  co{k\k2,  k 3};  if  not,  then  let  n  =  n  +  1  and  return  to  step  1;  if  yes,  then 
K  =  {k\k2,k3}  has  formed  a  simplex. 

Step  4:  Compute  the  sliding  vector  s  =  He(t)  and  o  =  [HB]~ls,  and  then  assign  the  simplex  U  =  {w1,  w2,  u3}  where  ul  = 
k1  •  | \s\ |,  to  be  the  control  vectors. 

Step  5:  Decide  to  which  open  region,  partitioned  by  the  simplex  control  vectors  U  =  {u\  u2,  u3},  the  o  belongs.  In  light 
of  the  definition  given  by  (7),  this  problem  can  be  solved  by  minimizing  the  following  objective  function: 


E,  = 


m+ 1 

-  E 

j=\J7i 


2juJ 


*  =  1,2, 


,m+  1, 


(20) 


subject  to  the  constraints  of  >  0.  The  MATLAB  command  Isqnonneg  can  be  used  to  resolve  this  plain  linear 
programming  problem  in  simulation.  If  there  exist  m  coefficients  >  0  such  that  Et  — >  0,  then  o  G  Qt. 

Step  6:  Apply  the  simplex  control  law  (8)  for  tracking  control  of  the  nominal  system  (14)  in  the  duration  of  the  sim¬ 
ulation. 

Step  7:  Compute  the  performance  index  rn+ 1  in  (18),  and  assign  the  optima  as  follows: 

If  (n  =  0)  or  (rn+ 1  <  r*)  then  T*  =  rn+ 1,  U*  =  Un+i  else  do  nothing; 
n  =  n  +  1,  and  then  return  to  step  1. 


Repeating  the  above  procedure,  in  finite  time,  one  can  get  the  optima  T*  and  U *,  which  minimize  the  performance 
objective  (18). 


4.  Simulation  study 

In  the  simulation,  the  objective  is  to  use  two  control  inputs,  u(t)  e  R2,  to  guide  the  state  x(t)  E  R 3  of  system  (1)  to 
match  a  prespecified  reference  signal  xr(t).  Note  that  this  is  much  more  difficult  than  the  tasks  achieved  before  [17,18]. 

The  reference  trajectory  xT(t)  is  specified  as  a  closed  orbit  corresponding  to  a  periodic  solution  of  the  unforced 
Chen’s  equation.  Let  the  parameters  of  system  (12)  be  a  =  45,  b  =  1.5  and  c  =  28,  so  that  system  (12)  generates  a 
periodic  solution  [15].  Starting  from  the  initial  state  jcr(0)  =  [—1.7570,  -1.9648, 7. 9743]T,  the  three-dimensional  phase 
portrait  and  the  time-domain  response  of  the  reference  jtr(f)  are  shown  in  Figs.  1  and  2,  respectively. 

For  the  purpose  of  controlling  the  nonlinear  uncertain  systems  (1)  by  sliding  mode  control,  the  desired  sliding 
manifold  was  firstly  constructed  via  the  method  proposed  in  [16]: 


Fig.  1.  The  reference  orbit  xr(t),  plotted  in  the  x3  —  jci  —  x2  space. 
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Fig.  2.  The  deterministic  time  series  of  the  reference  orbit  xr(t). 


S(e)  =  {e\ s(e,  t)  =  He(t)  =  0},  s(e)  E  R 2, 
where  e  =  x  —  xr  and 

[4.1680  7.1120  9.6191  ]  f  ) 

|_7.4080  9.7120  —3.4480  J '  (  } 

Next,  in  light  of  the  procedure  presented  in  Section  3,  the  optimal  simplex  control  vectors  were  sought  via  chaos 
optimization.  The  upper  bound  and  lower  bound  of  the  searched  variables  are  chosen  as  lbd7  =  —  50  and  ubdy  =  50,  and 
the  procedure  was  performed  for  about  300  times.  The  result  is  U*  =  {ul,u2,u3},  u!  =  k!  •  ||s||  and 

k1  =  [0.0778, 47.1  185]t, 

K2  =  [-49.992,  — 10.0773]t,  (22) 

k3  =  [26.2508,  — 3.5633]t. 

With  the  obtained  simplex  control  vectors  U* ,  the  tracking  control  performance  under  the  simplex  sliding  mode  control 
law  (8)  is  visualized  by  Figs.  3-6,  which  shows  that  the  trajectory  of  the  controlled  chaotic  system  is  steered  to  the 
reference  signal  with  satisfactory  performance. 


Fig.  3.  The  controlled  trajectory  x(t)  of  Chen’s  chaotic  system  to  the  reference  orbit  xT(t),  plotted  on  the  x3  —  jci  —  x2  space. 
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Fig.  6.  Tracking  performance  of  the  third  coordinate. 
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5.  Conclusions 

It  is  known  that  the  simplex  sliding  mode  control  has  made  several  important  improvements  over  the  conventional 
sliding  mode  control  for  robust  control  of  nonlinear  dynamical  systems,  and  yet  the  selection  of  the  magnitudes  and 
directions  of  the  control  vectors  in  this  method  is  usually  determined  by  trial-and-error  through  simulation  and/or 
experimental  testing.  In  this  paper,  a  systematic  controller  design  method  based  on  chaos  optimization  has  been 
developed  for  simplex  sliding  mode  tracking  control  for  affine-type  of  uncertain  nonlinear  systems  with  multi-inputs. 
This  method  provides  a  step-by-step  routine  for  constructing  the  simplex  control  vectors  for  the  nonlinear  sliding  mode 
control,  thereby  resolving  a  long-lasting  difficult  technical  problem  in  robust  control.  Computer  simulation  on  a 
nonlinear  uncertain  system  has  illustrated  the  effectiveness  of  the  proposed  design  method. 
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